CMB polarization power spectra contributions from a network of cosmic strings 
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We present the first calculation of the possible (local) cosmic string contribution to the cos- 
mic microwave background polarization spectra from simulations of a string network (rather than a 
stochastic collection of unconnected string segments). We use field theory simulations of the Abelian 
Higgs model to represent local U(l) strings, including their radiative decay and microphysics. Rel- 
ative to previous estimates, our calculations show a shift in power to larger angular scales, making 
the chance of a future cosmic string detection from the B-mode polarization slightly greater. We 
explore a future ground-based polarization detector, taking the CLOVER project as our example. In 
the null hypothesis (that cosmic strings make a zero contribution) we find that CLOVER should limit 
the string tension /i to Gfi < 0.12 x 10~ 6 (where G is the gravitational constant), above which it is 
likely that a detection would be possible. 



I. INTRODUCTION 

There is an increasingly strong observational case for 
the inflationary paradigm, including precise measure- 
ments of the cosmic microwave background (cmb) ra- 
diation. Many viable inflation models have been con- 
structed, under a wide range of physical theories, but a 
prediction in many interesting cases is that a network 
of cosmic strings [l|, |2( should exist after inflation. This 
includes all viable inflation models under grand unified 
theories (guts) below a certain complexity [3|l as well as 
brane inflation models in string theory • How- 

ever, the predicted cosmic strings would interact with ra- 
diation and matter via gravity and so source variations in 
the temperature of the CMB in addition to those seeded 
directly by inflation. If their energy per unit length is 
too large then they destroy the fit to data yielded by the 
inflationary component, but current measurements still 
allow the string component to the temperature power 
spectrum to be about 10% (dependent upon angular 
scale) @ i EH. Hence, while strings may be impor- 
tant for future cmb measurements, they produce at most 
a secondary component of the temperature anisotropies, 
which is shrouded by that of inflation and is therefore 
difficult to detect. 

However, the polarization of the CMB radiation con- 
tains information in addition to that from the tempera- 
ture variations and it has been shown that inflation con- 
tributes only weakly to an gula r polarization patterns of a 
so-called magnetic nature [ll| . CMB polarization is hence 
a particularly interesting means to detect anisotropy 
components that are sub-dominant in the temperature 
power spectrum, since they may still give the primary 
contribution to polarization patterns of this type. Two 
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important examples of such components are primordial 
gravitational waves and the case of primary interest here: 
cosmic strings [l2j • In this paper we present the CMB po- 
larization power spectra contributions from the latter, 
calculated for the first time using classical field theory 
simulations of a cosmic string network. The method for 
our calculations is discussed in more detail in [l3j . in 
which we limited the results presented to the temper- 
ature anisotropy. We reveal our polarization results in 
this separate publication because our discussion relies on 
results from [8| but also because the nature of previous 
temperature and polarization calculations are different 
and these are not only the first polarization results to 
stem from field simulations of cosmic strings, they are 
the first to use actual simulations of a cosmic string net- 
work. Finally we also wish to highlight that polarization 
offers exciting future observational possibilities and de- 
serves a detailed discussion. 



II. CMB CALCULATION METHOD 

A. Brief overview 

For this work we make use of the fact that the mea- 
sured CMB anisotropies are only of order 1 part in 10 5 
and therefore that the perturbations in the radiation 
and matter do not have a noticeable effect on the string 
dynamics. Simulations of cosmic strings can therefore 
be performed in a homogeneous Friedman-Robertson- 
Walker universe with negligible loss of accuracy. The 
energy-momentum tensor of the cosmic strings from these 
simulations is then used to source metric perturbations 
via the Einstein equations and so creates perturbations 
in the radiation and matter. The low level of the pertur- 
bations also implies that any interaction between these 
string-induced perturbations and the primordially seeded 
inflationary ones can be ignored, and therefore that the 
cosmic string contribution to the CMB can be calculated 
separately. With negligible coupling between them, the 
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two contributions are statistically independent and the 
total power spectrum is simply the sum of the power 
spectra from each of the two sources. 

The contribution from inflation has been well stud- 
ied and a number of codes exist which accurately solve 
the Einstein-Boltzmann equations and, using a modern 
desktop computer, results for both temperature and po- 
larization power spectra can be obtained in a few seconds. 
The cosmic string contribution is more challenging due to 
the non-linear string dynamics and the enormous range 
of scales involved in the problem. There is strong evi- 
dence that cosmic strin g ne tworks enter an attractor so- 
lution or scaling regime [H, Q QJ, QJ, Q3, Qj, [H M, HI 
in which the average string separation is comparable to 
and scales with the causal horizon size. However, the 
width of strings corresponds inversely to their energy 
scale, which may be ~ 10 16 GeV, and therefore the width 
is many orders of magnitude smaller than their separa- 
tion at times of importance for CMB calculations. The 
problem is hence not directly solvable with current or 
near future technology. Indeed, the dynamics of strings 
on scales much smaller than the horizon is not well un- 
derstood [H, H3, HH [HJ , with questions relating to the 
production of string loops and of small-scale structure on 
long strings. There is not even a consensus on whether 
gravitational radiation or high energy particle emission is 
the dominant means via which strings decay [3, H2| • 
These facts limit the reliability of gravitational radiation 
predictions but fortunately CMB anisotropies are sourced 
mainly by the distribution of long strings on scales close 
to the horizon, which are more easily simulated. 

The only previously published CMB polarization results 
for (local) cosmic strings @ (and the derived [HHHI) use 
the computationally rapid method of Albrecht et al. [24| , 
which is reliant upon two levels of approximation. Firstly, 
the Nambu-Goto approximation is applied, such that the 
strings are taken to have infinitesimal width. Secondly, 
rather than simulating a string network, the strings are 
represented by a series of straight and unconnected string 
segments, moving with stochastic velocities. The sub- 
horizon decay of strings is modeled by the random re- 
moval of segments such that the total string length varies 
as it would for a network according to the above scaling 
law. Then the angle-averaged energy-momentum tensor 
of the segments is used to source the cosmological per- 
turbations. The results for the temperature power spec- 
trum obtained from this unconnected segment model are 
in broad agreement with those obtained from actual sim- 
ulations of a connected string network, evolved using the 
Nambu-Goto equations of motion [25|. This gives con- 
fidence in the applicability of the computationally rapid 
method, however for the case of polarization, there are 
no published results from Nambu-Goto simulations with 
which to compare those found using the unconnected seg- 
ment model. 

However, the Nambu-Goto simulations themselves are 
not complete and questions have been raised over the 
accuracy with which they represent string loop produc- 




FIG. 1: A snapshot from an Abelian Higgs simulation in 
the matter era at a time when the horizon volume approx- 
imately fills the simulation box. The lines show the centers 
of the strings (found using the gauge-invariant phase- winding 
method of [271]) while the upper and lower surfaces highlight 
the additional presence of radiative decay, which must be in- 
cluded in an adhoc manner in Nambu-Goto simulations. The 
lower surface indicates regions of significant energy density 
due to the non- vacuum value of the Higgs field while the up- 
per surface shows regions of significant energy from the quasi- 
magnetic field in the model (see UU). Note, however, that the 
strings themselves make the primary contribution to both of 
these types of energy and the contrast is chosen to highlight 
the radiation contribution. For example, the circular pattern 
seen on the left in these slices is due to the recent collapse of 
a string loop just above the bottom of the simulation (and is 
seen in both slices due to the periodic boundary conditions). 



tion and the string decay [H [U, H H|. In [13] we 
presented the first calculations of the temperature power 
spectrum for (local) cosmic strings to stem from field- 
theoretic simulations. Those calculations included CMB 
polarization, as is required for the accurate calculation 
of the temperature anisotropy and we hence present our 
polarization results here with minimal discussion of the 
method employed. 

In short, by employing field simulations of the Abelian 
Higgs model, we have studied local U(l) strings with the 
string width resolved, accurately representing the asso- 
ciated microphysics, including a form of radiative decay, 
see Fig. [TJ While this means that the simulations were 
limited to very small length scales at very early times, the 
scaling property enables statistical results to be trans- 
lated to larger scales and later times and hence provide 
the information required for CMB calculations. The cor- 
responding extrapolation may be over many orders of 
magnitude, but it is justified by the behaviour of the 
strings in the simulations and it is not merely a theoret- 
ical assumption. 
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We now present greater depth of the calculation pro- 
cedure, in particular pointing out that scaling alone is 
not sufficient to enable current technology to perform 
the string simulations needed for CMB calculations, but 
readers who are familiar with [13j, or who are less con- 
cerned with these details, may wish to leave the following 
section for future reading. 

B. Greater detail 

Despite the use of the scaling property of cosmic 
strings, the simulations required for CMB calculations 
could not be performed directly. The width of an Abelian 
Higgs string is a fixed physical scale and hence lessens 
rapidly relative to the physical horizon size, which varies 
as r 2 in the radiation era and t 3 in the matter era, where 
r is the conformal time. Since scaling is broken at the 
radiation-matter transition, strings must be studied un- 
der both radiation and matter domination, but to do 
so over a sufficiently large range in r is very challenging. 
This is a particular problem because strings will not form 
in the simulations until the horizon is much greater than 
their width and even then, the strings will take some time 
to reach the scaling regime. 

Our solution to this problem was to raise the coupling 
constants in the Abelian Higgs action to time dependent 
variables such that the string width r grew slightly in 
physical (rather than comoving) units: 

rphys. oc a 1_s , (1) 

where a is the cosmic scale factor and s is a parameter 
that controls the growth. The true case has s = 1 while 
a fixed comoving width is obtained by setting s = 0. The 
later is perhaps the most straightforward to simulate and 
the authors of [19j had previously studied string scaling 
using equations similar to those obtained via our method 
with s — 0. The closest to the true case that proved possi- 
ble using 512 x 512 x 512 simulation lattices in the mat- 
ter era was s = 0.3 [7l| (with simulations performed using 
64 processors of the UK National Cosmology Supercom- 
puter |28J). However, the use of s + 1 1 causes a breach 
of the conservation law of the very energy-momentum 
tensor that sources the anisotropics, which may not have 
been an major problem for [l!| but is obviously a poten- 
tial problem for CMB calculations. Fortunately the effect 
can be investigated by simulating strings at a number of 
s values, including the true s = 1 case for the less diffi- 
cult radiation era. Accordingly, in [13j | we showed that 
s has only a minor effect on the string length behavior, 
the Fourier distribution of the energy-momentum ten- 
sor components and upon the CMB results. The use of 
s < 1 is therefore treated merely as a systematic source 
of uncertainty, comparable to the statistical uncertainties 
arising from realization-to-realization variations. Hence 
as a result of this procedure, the data required for accu- 
rate calculations of the CMB power spectra were obtained 
from the simulations, with the results below being taken 



from simulations with s = 0.3 in both radiation and mat- 
ter eras. 

If only power spectra are desired (rather than bi- 
spectra or maps), the unequal-time correlations of the 
energy momentum tensor T^: 

^(k,r,r') = (f AK (k ) r)f; i/ (k,r')) , (2) 

are all that is required to solve the Einstein-Boltzmann 

equations [UlUlIi While there are | 10 ( 10 + 1 ) = 55 
such correlation functions, each complex functions of 
3+1 + 1 = 5 input variables, they are heavily con- 
strained by scaling, statistical isotropy, causality and 
energy-momentum conservation. As a result they may 
be represented by merely 5 real functions of 2 variables: 
C a b{kT,kT') [Uj[2!J,[30]. For example, statistical isotropy 
implies that only the magnitude of k is important for 
any expectation value, while scaling implies that the spa- 
tial distribution, on average, scales with the horizon and 
hence only the product kr is important, rather than k 
and r individually [72[. Of the 5 scaling function, 3 
represent scalar degrees of freedom with ab equal to 11, 
12, 22 (with the 21 case given by r <-> r'), while vec- 
tor and tensor degrees of freedom give one function each. 
These unequal-time scaling functions (calculated under 
both matter and radiation domination) encode the infor- 
mation required to accurately calculate the power spectra 
contributions as a result of the sourcing T M „ . 

However, in order to solve the Einstein-Boltzmann 
equations it is convenient to re-express these scaling func- 
tions, which like U\^ u .u are quadratic in T^ v , as linear 
functions c„(fcr) |13.l3ll|: 

C(kT, kr') = J2 Kcn(kr) c n (kr'). (3) 

n 

In the numerical case, C is known only at discrete val- 
ues of kr and kr' and the above equation corresponds 
to a decomposition of the real, symmetric matrix C into 
its eigenvectors c„ and eigenvalues A n . The eigenvectors 
then represent linear combinations olT^ Vl which source 
metric perturbations via the Einstein equations^ These 
were fed into a modified version of CMBEASY (32j, and 
the CMB power spectrum contribution from a particular 
eigenvector was found, with the total string contribution 
being the sum over these. Modifications to CMBEASY in- 
cluded the additional presence of the source terms and 
also cosmological perturbations of a vector nature, which 
decay in the standard inflation case but are continuously 
sourced by cosmic strings. These were evolved using 
gauge-invariant perturbation equations obtained via the 
Hu and White total angular momentum method [33j . 

III. RESULTS 

The results for both temperature and polarization 
power spectra are shown in Fig. [51 and compared to 
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FIG. 2: The CMBtemperature and polarization power spectra 
contributions from cosmic strings (black), inflationary scalar 
modes (gray, solid) and inflationary tensor modes (gray, dot- 
dashed). For the case of the te cross correlation, positive 
correlations are shown as solid lines and anti-correlations are 
shown as dashed lines, except that is for the inflationary ten- 
sor component for which the sign is not indicated here. 



the inflationary contribution, calculated here using CAMB 
[34| . The normalizations of the three quantities shown 
in the figure: the string contribution, the inflationary 
scalar contribution, and the inflationary tensor mode; are 
each free parameters. The inflationary tensor (or grav- 
ity wave) contribution has yet to be discovered, with the 
scalar mode alone giving an excellent fit to current data. 



Hence in the four plots shown, the inflationary scalar 
component is shown with its normalization A 2 and other 
cosmological par ameters chosen to match current CMB 
data dfl l3ol l37l [H, H^], without including strings or 
inflationary tensor modes (see endnote [73j for parame- 
ters). 

The normalization of the string contribution is set at 
the 95% upper bound allowed by current CMB data, as 
calculated by the current authors in Q , using the above 
CMB data together with parameter estimates from non- 
CMB experiments: h = 0.72 ± 0.08 and f2 b /i 2 = 
0.0214 ± 0.0020 HH. This bound assumes negligible pri- 
mordial tensor modes (although it is not especially sen- 
sitive to them) and corresponds to an 11% contribution 
to the temperature (tt) power spectrum at I = 10 or 
fractional contribution /10 = 0.11. For this string com- 
ponent, the normalization is related to the energy-scale 
of the associated theory, which determines the string ten- 
sion /z. It is conventional to express this in terms of G[i 
(where G is the gravitational constant) with the power 
spectra proportional to (Gfi) 2 . The results plotted here 
correspond to a value Gti = 0.7 x 10~ 6 Q (and also 
h = 0.72, fl h h 2 = 0.0214, n A = 0.75 0, spatial flatness 
and t = 0.1). 

In a similar manner, the inflationary tensor mode nor- 
malization is at the 95% upper bound set by the same 
data and method that we employed in Q to constrain 
the string component, without including strings while 
calculating the tensor case. This corresponds to an in- 
flationary tensor contribution to the temperature (tt) 
power spectrum at multipole £ — 10 of 15% or a ratio of 
tensor-to-scalar primordial perturbation power spectra of 
r = A 2 /A 2 = 0.36 (at comoving wavevector fco = 0.01 
Mpc" 1 ). In this determination we used the single-field 
consistency equation [43| to give the tensor spectral tilt: 
nt = —r/8, which implies n t = —0.045 for the case plot- 
ted. 

When interpreting the results, it is useful to consider 
the mechanism involved in creating CMB polarization: 
Thomson scattering in the presence of a quadrupole in- 
tensity anisotropy, as seen by the scattering particle. 
This results in differing contributions being given to per- 
pendicular polarization directions and therefore partially 
polarized emission. However, in the presence of very 
strong Thomson scattering, a quadrupole moment can- 
not be set up and hence there is no contribution until the 
universe begins to deionize. But as this occurs the fre- 
quency of scattering decreases, soon becoming negligible. 
Hence CMB polarization stems primarily from a short pe- 
riod around deionization and therefore the polarization 
spectra are dominated by small angular scales. However, 
the radiation from stars results in partial reionization for 
more recent times, giving a small peak at large angular 
scales. These two features are clearly seen in all polar- 
ization graphs, for both strings and inflation, although 
the inflationary tensor modes do not contribute to very 
small scales. 
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A. TT spectrum 

Before discussing the individual polarization spectra, 
the tt power spectrum component from strings deserves 
a brief comment here. The string contribution is a broad 
peak between I w 40 and £ « 600 (75% of peak value), 
without the acoustic oscilations that are seen in the data 
and the matching inflationary scalar mode. In the in- 
flationary case there are frozen super-horizon perturba- 
tions, which for a particular length scale enter the grow- 
ing horizon and begin to oscillate at the same time. This 
gives a high degree of temporal coherence. In the string 
case, perturbations are continuously sourced within the 
horizon by objects which themselves experience complex, 
non-linear dynamics. As a result no oscillations are seen 
in the plot for strings. 



B. EE spectrum 

Unlike the decomposition of a spin-1 vector field into 
curl-free and divergence-free parts, the decomposition of 
the spin-2 polarization field into so-called Eand Bmodcs 
is based upon the second angular derivatives of this 2D 
field. We first address the E polarization power spec- 
trum, which measures polarization patterns in which a 
minimum in the polarized intensity is accompanied by 
the greatest gradient being in a direction parallel or per- 
pendicular to the plane of polarization. Only then will 
we discuss the B mode, for which the greatest gradient is 
at 45° to the polarization plane. 

The cosmic string component to the EE power spec- 
trum is, very approximately, of the same form as that 
from inflationary scalar modes, although with much 
smaller oscillations at small scales. Since strings source 
perturbations most strongly on scales that are well within 
the horizon, whereas inflation seeds super-horizon per- 
turbations that oscillate after they enter the horizon, the 
smaller oscillations are not surprising — nor is the ad- 
ditional shift to higher multipoles. Some oscillations are 
expected since the electric polarization is heavily domi- 
nated at small scales by scalar modes (which oscillate), 
but decoherence almost completely erases the oscillations 
from the EE power spectrum. 

The logarithmic vertical axes of these four plots share a 
common log-interval per unit length and therefore a given 
distance traversed upwards corresponds to an increase by 
the same factor on each plot. It is hence clear that the 
amplitude of the string EE component relative to that 
from inflationary scalar modes is comparable to that of 
the temperature power spectrum, but since temperature 
data is likely to be significantly more precise than polar- 
ization data for the foreseeable future, measurements of 
the EE data is not likely to place direct constraints upon 
cosmic strings for some time. 



C. BB spectrum 

The situation is quite different, however, for the 
BBpolarization power spectrum. Inflationary models only 
create scalar and tensor modes, while cosmic strings also 
create relatively large vector modes. The primordial 
scalar mode contributes to BB polarization mainly (74| 
through the partial conversion of the EE spectrum to a 
BB contribution, which is due to gravitational lensing 
by matter perturbations [75| . Therefore the contribu- 
tion from cosmic strings dominates the BB spectrum for 
i = 150 — 1000, despite them being constrained to be 
sub-dominant in the temperature power spectrum. At 
the current upper limit on the string normalization, the 
bb contribution is roughly 3 times larger than the in- 
flationary contribution in the range £ = 150 — 500 and 
despite it being likely that any real string contribution 
may have a lower normalization, a quite large reduction 
is required to stop them from dominating these scales. 

At low multipoles the reionization peak yields an even 
greater dominance over the lensing signal, however, the 
inflationary tensor modes may contribute significantly to 
these scales. Indeed, if strings contribute at the 95% up- 
per bound quoted above and make up 11% of the temper- 
ature power spectrum at £ — 10, while inflationary tensor 
modes contribute 15% (although strictly, these two upper 
bounds should not be employed simultaneously), then 
the tensor component would be between 2 and 8 times 
larger over the range 20 < £ < 100, but about equal to 
that from strings for £ « 10. As a result, if strings con- 
tributed at the level plotted but r < 0.05 (as would be 
the case for supersymmetric hybrid inflation models (lo| ) 
then strings would dominate at all £ < 1000. 

D. TE spectrum 

The te cross-correlation spectrum (between the tem- 
perature and electric polarization mode) shows signifi- 
cant oscillations in the string case, although decoher- 
ence suppresses them relative to the inflationary case. 
Both the inflation and string contributions are constant 
to within an order of magnitude for large scales and they 
both oscillate between anti-correlations and positive cor- 
relations for 100 < £ < 1000. Note that the sign conven- 
tion we chose for the te spectrum is that by the wmap 
team. 

As is clear from the figure, the te spectrum contri- 
bution from cosmic strings is generally less significant 
relative to the inflationary contribution than in the tt 
or bb spectra. Although the oscillations seen for cos- 
mic strings are 180° out of phase with the those from 
the inflationary scalar mode, which may assist the use 
of te data in constraining cosmic strings, it is evident 
that the te spectrum not especially important for con- 
straining the cosmic string component directly. Future 
te data will however, further constrain the inflationary 
scalar contribution to all spectra and reduce the cosmo- 
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logical parameter degeneracies which allows it to change 
and accommodate large cosmic string (or other topologi- 
cal defect) components while still fitting the data [3,[4J]. 

E. Uncertainty estimation 

Finally, it is of course the case that all these results 
incorporate, at some level, uncertainties relating to the 
random initial conditions used in the simulations, finite 
volume effects and the use of s < 1 (see Sec. HI|) . Our 
investigations into these and other possible uncertainties 
are described in the appendix. 



IV. COMPARISON WITH PREVIOUS LOCAL 
STRING RESULTS 

The only previously published calculations of the po- 
larization power spectra contributions from local cosmic 
strings 0] are based on the unconnected segment model. 
This model has a number of free parameters which were 
fixed by reference to Nambu-Goto simulations [ItJ and 
another phenomenological model, the velocity-dependent 
one- scale model [II EI 111 . In we pointed out that 
our simulations lead to temperature power spectrum re- 
sults that are more biased to large scales than those of the 
unconnected segment model |47j . with the vector com- 
ponent being of greater importance and having a quite 
different form. Reducing the so-called wiggliness param- 
eter in the unconnected segment model can to a certain 
extent boost the vector contribution relative to the scalar 
component, as can be seen in Fig. 1 of 0, but the model 
cannot include velocity correlations (or the angular mo- 
mentum of the decay products) which source the vector 
mode. 

In this discussion of our polarization results, we addi- 
tionally note that there is a similar shift in angular scale 
in the vector-dominated bb spectrum. Our string simu- 
lations give a bb peak at £ = 600, with half of the peak 
value given at I = 300 and 1100. On the other hand, the 
simpler model gives a peak at I = 800 with the half of 
the peak value given at I = 400 and 1500. This may be 
important observationally because our results highlight 
a more significant difference in form between the cosmic 
string contribution and that from gravitational lensing. 

The amplitude of the power spectrum also differs be- 
tween the two calculation methods, with the unconnected 
segment model giving values for both peaks that are ap- 
proximately 10 times greater than ours for a given Gp, 
value. However, the normalization of the temperature 
power spectrum is also larger for that model and given 
constraints provided by current temperature data this ef- 
fect is largely removed, with the factor of 10 becoming a 
mere doubling. 

Hence it would appear that the computationally rapid 
unconnected segment method does yield qualitatively 
correct polarization results but with quantitative differ- 
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FIG. 3: A comparison of the EE (dashed) and bb (solid) 
power spectra contributions from local cosmic strings (black) 
and global textures (gray). Each are calculated at h = 0.72, 
Q b h 2 = 0.0214, Q A = 0.75 and r = 0.1, with the normaliza- 
tions set to yield fractional contributions to the temperature 
power spectrum of fio =0.11 for strings and /io = 0.13 for 
textures, which correspond to the 95% upper bounds from 
CMB data. 

ences relative to the results from our more complete sim- 
ulations. 



V. COMPARISON WITH GLOBAL DEFECTS 

Global defects result from the spontaneous breaking of 
global symmetries in the early universe, and share the 
important scaling property with cosmic strings. Unlike 
their local counterparts, global defects do not localize 
their energy into the defect cores and, for example, in the 
case of global strings, the strings themselves do not need 
to be resolved in the simulations. Hence, without the 
need for our string width control formalism (see Sec. HT| . 
temperature power spectra components from field the- 
ory simulations have previously been published for global 
texture s filj l29l liil . , global monopoles [29| and global 
strings [29|]. Polarization spectra for these three global 
defect types have also been presented previously [49( and 
we have applied our unequal-time correlation software 
and modified CMBEASY code to the case of global tex- 
tures to provide one means of checking our results against 
independent polarization calculations, and find excellent 
agreement. 

This also enables us to consider the differences in the 
polarization signals from local strings and global tex- 
tures, both calculated at modern cosmological parame- 
ter values and shown in Fig. [3l In the figure we have 
normalized the global texture contribution to match the 
95% upper bound found in [44[ using CMB data and the 
above measurements of h and fl^h 2 . For the tempera- 
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ture power spectrum (see [13|), the key difference is that 
the textures yield a large bias toward low multipolcs in 
comparison to local strings, giving a broad peak between 
£ = 20 and £ = 300 (75% of peak value). This trend can 
also be seen in both the EE and bb polarization spectra, 
with the later peaking at £ = 200 in the texture case, 
rather than at £ = 600 as it does for strings. 

This can be understood as follows. The decay of topo- 
logical defects and the smoothing of the associated fields 
is limited by the causal horizon (or more precisely the 
distance that light could have traversed since inflation 
ended) . In the global case there are large-scale gradients 
which constitute a form of potential energy. The dynam- 
ics of the system therefore rapidly removes the gradients, 
smoothing out the textures on scales just a little smaller 
the horizon and textures only contribute to energy and 
momenta for kr < 10 (13 . l4a ]. In the local string case 
however, a gauge field is present which can accommo- 
date the large-scale gradients such that they are effec- 
tively smoothed and contribute little energy. The strings 
themselves represent spatially extended potential energy 
contributions and the Hubble-damped dynamics tends to 
reduce the string length, but this process is less efficient 
so strings persist within the horizon. Hence, in contrast 
to the global case the Abelian Higgs fields contribute en- 
ergy and momenta for all fcr < 100 (l3| . Hence, the CMB 
anisotropics sourced by the energy-momentum tensor of 
the local string model are more biased to smaller scales 
than for global textures, and global defects in general. 



VI. CURRENT AND FUTURE CMB 
MEASUREMENTS 

These first polarization results from string network 
simulations are very encouraging since although a finite 
B-mode polarization has yet to be detected, there are a 
number projects in preparation that aim to measure B- 
mode polarization with high precision. Current CMB po- 
larization data gives merely an upper bound on the bb 
spectrum, with for example the full-sky wmap project 
providing a upper limit of 3 x 10~ 2 nK 2 derived from 
2 < £ < 12. This is only a factor of 3 above the plotted 
tensor mode prediction for r = 0.36, but about 30 times 
the mean value for G/i = 0.7 x 10~ 6 cosmic strings over 
this range. On smaller angular scales, the ground-based 
DASI and balloon-borne boomerang detectors place lim- 
its on the bb spectrum of a few /j>K 2 for multipolcs 
200 < £ < 600 [H [EH while the ground-based CBI 
project similarly constrains multipoles 400 < £ < 1700 
(52 |. These experiments were designed to detect to the 
EE and te spectra, and were successful in this aim, but 
the detection of B-mode polarization will require the next 
generation of sub-orbital CMB experiments and, on the 
largest angular scales, the forthcoming Planck satellite 
mission. 

Planck is scheduled for launch in 2008 and while its 
polarization sensitivity and angular resolution are signif- 
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FIG. 4: Estimated uncertainties from the future CLOVER 
project under the null hypothesis, that is that the true un- 
derlying model contains neither cosmic strings nor primordial 
tensor modes. The thick solid line indicates the underlying 
model, while the squares show the ^-binning of the data and 
the lcr uncertainties estimated. The thin solid line shows the 
total power spectrum from an /io = 0.0035 model while the 
dashed thin shows the string component in that case. 



icant improvements over WMAP, it is not specifically de- 
signed to detect very small B-mode signals. It will, how- 
ever, provide a measurement of the bb spectrum at large 
angular scales that will not be superseded until a subse- 
quent CMB satellite is launched, but this is more relevant 
for the detection of the primordial tensor mode than it is 
for cosmic strings. On smaller scales, Planck will struggle 
to detect the gravitational lensing signal, with very large 
uncertainties even after binning across large ranges in £. 
However, if strings contribute at the level plotted, they 
will be detected by Planck's BB measurements, as what 
might be otherwise interpreted as a gravitational lensing 
signal that is 3 times larger than expected. Of course, 
they would probably also be significantly detected as a 
result of the precise TT, te and EE data that Planck will 
provide. With full-sky tt data that is cosmic variance 
limited to multipoles possibly as high as £ = 2500 and the 
EE power spectrum measured with ~ 10% accuracy out 
to £ > 1000, the cosmological parameters will be heavily 
constrained. Then, the tt data should be sensitive to 
smaller string contributions than the current data. 

However, planned sub-orbital projects have the po- 
tential to be more sensitive to the B-mode polarization 
than the Planck satellite and, while they will not pro- 
vide data for the very largest angular scales, it is the 
range 100 < £ < 700 that is most interesting for cos- 
mic strings. These projects can deploy new technology 
more rapidly than a satclite mission and they can also 
choose the view fields to minimize foreground contami- 
nation. With the list of such projects including [53j, |54| : 
ebex (first flight 2008), PolarBeaR (operational in 2008), 
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QUlET(deployment in early 2008), Clover (operational 
in 2009), PolarBeaR II (operational in 2010), and spider 
(first long-duration flight 2010); there should be a wealth 
of B-mode data available in the relatively near future. 

We will consider just one of these as an example, 
Clover: a ground-based project that will measure the 
B-mode polarization over the range 20 < I < 1000 and is 
scheduled to become operational in 2009. The aim is to 
eliminate detector noise and foreground contamination 
to yield the precision required to detect primordial ten- 
sor modes even if r is as low as 0.01 j54|. However, it 
will also characterize the high I range and will be sensi- 
tive to string components that are very small in the tt 
spectrum. 

Given the sensitivity, resolution, and sky coverage of 
a CMB project, the likely uncertainties in the measured 
Cf B are approximately given by [11, 55, 56]: 



(21 + l)/ sky A£ 



(Ct + my 



(4) 



The first term is due to (Gaussian) cosmic variance, while 
the second is the thermal noise: 



Nt — #FWHM "pixel ex P 



1(1 + 1)#FWHM 



8 In 2 



(5) 



The C£over project is expected to cover a fraction 
/sky = 0.025 of the sky which means that measurements 
at each £ are not independent, with correlations extend- 
ing roughly A£ = 40. The equation assumes a Gaussian 
beam profile, which we take to have #fwhm = 8' and 
then gives map pixels representing solid angle ss $fwhm- 
We assume successful foreground removal and take the 
thermal pixel noise to be o~ p i xe i = 0.65 fxK, which cor- 
responds roughly to that expected for Clover after 2 
years of observation. We then suppose that both pri- 
mordial tensor modes and cosmic strings make a zero (or 
negligible) contribution — the null hypothesis — and in- 
sert the corresponding bb power spectrum into Eqn. [H 
with the results shown in Fig. 2] 

To make a crude estimate of the likely upper bound 
on the string contribution under the null hypothesis, we 
perform a one-dimensional likelihood analysis, varying 
/io- We reduce the normalization A% of the primordial 
scalar component, such that it is proportional to (1 — ,/io) 
but keep all other cosmological parameters fixed. This 
is partially justifiable since, by the time of the CLOVER 
two-year data release, the tt, te and EE data from the 
CMB, combined with other independent data, will im- 
pose tight limits on them under the standard empiri- 
cal models. Also the dependence of upon /io, which 
strictly could deviate from the chosen form for a fit over a 
range of multipoles, was found in [8| to be approximately 
(1 — /io) for current data. We then find that the 95% up- 
per bound on the string contribution stemming from the 
BB data will be f w < 0.0035 ± 0.0012, with the central 
value being a factor of 30 less than the limits given by 
current data Q. Note that there is a statistical uncer- 
tainly in this result because the measured data will be 



spread around the underlying model according to the er- 
rors indicated in the figure, and different realizations will 
yield different upper bounds. The corresponding string 
tension limit is G/x < (0.12 ± 0.02) x 10~ 6 . 

This is merely a first step towards a full investigation 
of possible future data, which would involve many dif- 
ferent projects and a full multi-parameter analysis, but 
the approximate result here suggests that projects such 
as CLOVER will be able to either detect cosmic strings 
from their B-mode contribution or place stringent limits 
on cosmic string scenarios. 

Note that in contrast with [T2], we have not inves- 
tigated the cleaning of the gravitational lensing signal 
from the CLOVER measurements, which in the standard 
case is partially achievable via the non-Gaussianity in- 
troduced by the lensing [H?], [H, This is because 
the cosmic strings will introduce non-Gaussian signatures 
themselves and the non-Gaussianity from strings is not 
well characterized at present. It should also be noted 
that our results are not greatly changed by the inclusion 
of a small tensor contribution, since tensor modes con- 
tribute only to the largest scales, while supersymmetric 
hybrid inflation models [10], for example, do predict the 
negligibly small tensor mode assumed here. 



VII. CONCLUSION 

We have presented the first CMB polarization calcu- 
lations to involve simulations of a (local) cosmic string 
network, rather than using unconnected string segments 
with stochastic velocities. We have demonstrated, both 
here and in a sister paper 13] in which we presented 
the method and temperature results, that field-theory 
simulations of local cosmic strings can be employed for 
this purpose with current computational facilities. Our 
results show that the computationally simpler uncon- 
nected segment method, which has been commonly used 
by other authors, gives results that are qualitively accu- 
rate. This had previously been confirmed only for the 
temperature power spectrum. However, our results do 
show a greater bias to large angular scales. 

Importantly, we confirm the prediction of a large string 
contribution to the B-mode polarization power spec- 
trum, even for small contributions to the temperature 
anisotropics, with current upper bounds from the tem- 
perature data allowing a possible dominance of the B- 
mode data by cosmic strings. Through simulating future 
data for the CLOVER project we have shown that the 
likely B-mode measurements appear sensitive to cosmic 
strings with fractional contributions to the temperature 
power spectrum at I = 10 as low as /io < 0.0035, which is 
a factor of 30 tighter than the current CMB bounds. This 
corresponds to an upper bound on the string tension of 
Gfi < 0.12 x 10~ 6 . 

This is very encouraging because a detection of cosmic 
strings would not only be interesting in itself, it would 
open a powerful window upon early universe physics. 
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We have included here a comparison between the pre- 
dictions for local U(l) cosmic strings and global defects, 
but different string models would of course also yield dif- 
ferent results. For example, it would be interesting to 
perform CMB calculations for semilocal string networks 
IsElHjmifiiHliSISfl, and also for simple models that 
yield Y-junctions between strings, such as the U(l) xU(l) 
dual Abelian Higgs model [67], 168J] or the non-Abelian 
SU(2)/Z 3 model of In string theory, composites of 
D- and F-strings are possible with Y-junctions where 
they separate into their constituents. If differences be- 
tween the B-mode spectra of superstring-inspired models 
and traditional U(l) strings could be established, then 
the detection of one or other type of cosmic string might 
represent an exciting observational test for string theory. 

Finally, we note that the CMB power spectra are not 
the only means to constrain or detect cosmic strings. 
Strings are likely to produce significant non-Gaussianity 
in the CMB and are not fully described by power spectra 
alone. In fact, cosmic strings would create discontinu- 
ities in CMB maps of size AT/T ~ 13G/x fl. These will 
only become directly noticeable once the pixel noise is 
of the same order as the discontinuities and the beam 
width of the microwave detector is sufficiently small for 
the inflation-induced CMB temperature variations on sub- 
pixel scales (which contribute an effective noise) to be 
similarly reduced. However, a statistical approach may 
be more sensitive to the string non-Gaussianity since it 
can draw information from a large number of pixels and 
hence is an interesting avenue for future research. 
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APPENDIX A: UNCERTAINTY ESTIMATION 

In this appendix we detail estimates of the uncertain- 
ties in our results, however, we also refer the reader to 
[l3| for a more in-depth discussion of both the method for 
calculating the power spectra and the limitations therein. 

An obvious source of uncertainty in our results stems 
from the fact that the initial conditions for the string 
simulations are randomly generated. Therefore the sta- 
tistical results taken from the simulations and used in 
the CMBEASY calculations will differ from realization to 
realization. We averaged the simulation results over 5 
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FIG. 5: The effect of s on the results compared to the esti- 
mated uncertainties from realization-to-realization variations, 
with shaded areas showing the 1 and 2-cr regions at each mul- 
tipole. Note that correlations extend over large multipole 
ranges, that there are statistical uncertainties for each s, and 
that T is the mean CMB temperature. 



realizations (in each of the matter and radiation eras), 
which amounted to 625 CPU-days of processing time on 
the UK National Cosmology Supercomputer [2g|. Hence 
we may estimate the variance in our power spectra by 
additionally applying the modified CMBEASY code upon 
the results from individual realizations, giving the results 
shown in Fig. \5\ In the case of the TE power spectrum, 
the estimated standard deviation in the mean over the 5 
realizations is about 20% for I > 1000 and is less than 
10% for £ < 100 but with larger relative values in the in- 
termediate regime where zero-crossings are common. In 
the EE spectrum, we find relative values between 4 and 
11% over the range 2 < £ < 3000 while in the BB spec- 
trum we find between 5 and 10% for £ < 1000, before this 
value increases roughly linearly with I to around 30% at 
£ = 3000. 

As discussed in Sec. [Hi the computational limitations 
and the rapid growth of the horizon size relative to the 
string width meant that the true Abelian Higgs dynamics 
had to be approximated such that the string width grew 
in physical units in proportion to a 1_s . The closest to the 
true s = 1 case that was achievable in the matter era was 
s = 0.3 and hence we investigated the effect of s upon our 
results in order to estimate the corresponding systematic 
error. It can be seen in Fig. that in the TE case there is 
a large change from the s = 0.3 and s = 0.2 cases and yet 
the s — case lies very close to that for s = 0.3. Further 
it should be noted that there are statistical uncertainties 
on all three lines and the results are consistent with the 
seen deviations between the three cases being largely due 
to the natural realization-to-realization variations rather 
than to any sensitive s dependence. 

More evidence in support of this conclusion comes from 
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the behavior of the strings in the radiation era, which can 
be studied under s — 1 for a meaningful period of time. 
The results for both the string length density and the 
equal-time correlation functions of the energy momen- 
tum tensor are shown in [13l ] to deviate from the s = 
case by an amount of the the same order as the statisti- 
cal uncertainties, with non-resolvable differences between 
the s — 0, s = 0.2 and s = 0.3 cases. Hence any system- 
atic errors in the te power spectrum results due to the 
use of s < 1 are likely to be of the same order as the 
measured statistical uncertainties, and this is also true 
for the EE and bb spectra. 

A further source of uncertainty stems from the possi- 
bility that the approximate scaling behavior seen is not 
actually the true behavior in the large r limit. While 
in principle this could be because the string decay prod- 
ucts included in the simulations are those inherent to the 
Abelian Higgs model rather than gravitational radiation 
which may in fact be more important. Unfortunately a 
detailed knowledge of the string dynamics on scales much 
smaller than the horizon is required before the nature of 
the primary decay products of cosmic strings are under- 
stood and there is no numerical estimate that is readily 
achievable for this factor. However, even ignoring the 
possibility of a important back reaction from gravita- 
tional radiation, the seen scaling behavior may still be 
in error. 

Under the expected scaling behavior, the mean (co- 
moving) length of string in per (comoving) horizon vol- 
ume Vh should be proportional to r, while Vh is propor- 
tional to r 3 and hence: 

£=ygoc.. (Al) 

However, this was not the scaling behavior seen and in- 
stead, £ was found to vary as: 

£oc(T-7i=o), (A2) 

in the scaling era [HI . While this is not a problem, since 
at the late times required for CMB calculations the offset 
T£=o is negligible, it does mean that every r used in the 
calculation of the unequal-time scaling functions must be 
replaced by r— r^ = o m order for them to be representative 
of the late time behavior. If there is any curvature in the 
supposedly linear region, then the measured t^—q is be a 
function of the r range over which a linear fit is applied. 
For the final results presented herein we fit over the range 
64(/)^" 1 < r < 1280Q 1 , where (ftp is the energy scale in the 
Abelian Higgs model (see |13j). By varying the upper 
limit of the fitting range a ball-park estimate of the effect 
may be obtained, with the results shown in Table U for 
I = 300, which is roughly where the contribution to the 
BB spectrum from strings is most significant. 

An additional systematic uncertainty arises from the 
limited range of time ratios at which the unequal time 
correlation functions can be calculated. As explained in 
[lH], the data is taken from the simulations over the pe- 
riod in conformal time: 64</>o < r < 160aV This is the 
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-44 


-5.5 
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Matrix size M = 512 — » 256 


+8.5 


+0.15 


-1.4 
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+25 


-13 


-17 



TABLE I: Investigations into the uncertainties in the CMB 
power spectra contributions from cosmic strings at £ = 300. 
The estimated uncertainty, or the responses to changes in the 
calculation procedure, are shown as percentages. 



time period after the strings have formed and are scal- 
ing to well within the statistical uncertainties but before 
the simulation becomes noticeably aware of the periodic 
boundary conditions. It gives a maximum time ratio of 
160/64 = 2.5, however, the offset scaling law and the in- 
corporation of Tf=o decreases this value to around 1.8. 
Hence, we cannot study the scaling regime over large 
ranges of time. Fortunately, the unequal-time correla- 
tion functions decay for large and small time ratios (see 
[H[) and hence the most important data can be taken 
from the simulations that were possible and we took the 
correlation functions to be zero for more extreme time 
ratios than we were able to study. 

The error instilled by this can be investigated by fur- 
ther zeroing known regions of the correlation functions, 
above a time ratio i? max and below a time ratio R^ x - 
We explored the change in the power spectra for trun- 
cations with R max = 1.6, 1.5 and 1.4, which led to the 
changes shown in the table. The te results are very sensi- 
tive to i?max and a further investigation may be required 
before a comparison of these results against, for example, 
Planck data can be reliably performed. A similar situa- 
tion exists for the EE data over certain £ ranges: £ ~ 20 
and £ > 1000, however, again it will be some time before 
data exists that will necessitate te and EE calculations 
for strings of greater precision. Fortunately, the more 
important bb spectra are robust, with the present re- 
sults appearing more than adequate for comparisons to 
BB data from, for example, CLOVER. 

A further source of uncertainty arises as a result of 
possible numerical errors during the eigenvector decom- 
position phase. The simulations output data for a partic- 
ular timestep r against a reference time r' for a series of 
k values and hence output C(kr,T /t') at discrete values 
of its two inputs. However the eigenvector decomposi- 
tion requires the data to be available at discrete kr and 
kr' to form a matrix of size M x M. There is hence an 
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incorporation of numerical errors during the required in- 
terpolation procedure that is more important for smaller 
M. Further, the resulting eigenvectors have a kr reso- 
lution that is dependent upon M (if linear kr spacing is 
used over a given kr range) and interpolation is carried 
out to provide data for the modified CMBEASY code at a 
given timestep r when solving a given k mode. We hence 
explored a reduction in M from the value of 512 that is 
used for our primary results, giving the changes shown 
in the table at I = 300 and changes less than about 10% 
in the EE and BB spectra for all t < 1000, and in the TE 
spectra for £ < 100 (before the oscillating regime begins). 
Note that larger values of M necessitate a larger number 



of eigencontribution terms to be included inorder to yeild 
convergence of the power spectra sums to well below the 
other error estimates. 

Finally, we investigate the possible errors in our mod- 
eling of the radiation-matter transition, at which scaling 
is broken and we interpolate between in the ei genv ectors 
(and eigenvalues) calculated in each era (see [lj]). An 
overestimate of the effect can be made by replacing the 
radiation era data with that from the matter era, yield- 
ing the results shown in the table for I = 300. However, 
the change is much less important at large angular scales 
where the impact of the radiation era is smaller. 
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